OPERATORS ASSOCIATED WITH A STOCHASTIC DIFFERENTIAL EQUATION 
DRIVEN BY FRACTIONAL BROWNIAN MOTIONS 



FABRICE BAUDOIN, LAURE COUTIN 

Laboratoire de Probabilites et Statistiques 
Universite Paul Sabatier 
31062 TOULOUSE Cedex 9 France 
fbaudoin@cict.fr, coutin@cict.fr 



Abstract. In this paper, by using a Taylor development type formula, we show how it is possible to 
associate differential operators with stochastic differential equations driven by a fractional Brownian 
motion. As an application, we deduce that invariant measures for such SDEs must satisfy an infinite 
dimensional system of partial differential equations. 



Contents 



4. Expectation of iterated integrals of the fractional Brownian motior 



1. Introduction and main result 



^^^^^^xr^oMc^ewlorjrn^n^^i^m^^^m^s^f^^ 



^^^4r^icati(^i^^^c^t^f^^f^Tr^mMr^jri£asrrrei 



6^^^Arjp^i^^^^0nj^uj^Tj3a^T^icon; 



6.1. Basic definitions and propertie; 
^^^^^^^^^^^^^^^j3ath|_ 



J 



6 i 3^ i ^rrtejjy^ti£r^i i hcOTv_d 



6.4. Ito maps : rough path with 2 < v < 3 



References 



1 

4 
6 
9 
17 
17 
17 
19 
22 
23 
26 



1. Introduction and main result 
A ci-dimensional fractional Brownian motion with Hurst parameter H G (0, 1) is a Gaussian process 

B t = {Bl,...,Bf), t>0, 
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where B 1 , B d are d independent centered Gaussian processes with covariance function 

R&a) = \{* B + * B -\t-a\* B ). 

It can be shown that such a process admits a continuous version whose paths have p finite variation for 
1/p < H. Let us observe that for H — i, B is a standard Brownian motion. 

In this paper, we are interested in the study in small times of stochastic differential equations on W 1 
(1.1) X?=x fi +Y,fvi{X?)dB\ 

where the V^s are C°°-bounded vector fields on R n and B is a d dimensional fractional Brownian motion 
with Hurst parameter H > -|. 

Let us recall that a smooth vector field V on R™ is simply a smooth map 

V : E™ -> M" 

a; i-> (^i(x), ...,v n (a:)). 

It defines a differential operator acting on the smooth functions / : R™ — * K as follows: 

(F/)0r)=f>(:r)£. 

i— 1 

With this notation, we observe that V is a derivation, that is a map on C°°(R n , R), linear over R, satisfying 
for f,g € C°°(R",R), 

If if > |, the integrals 

[ Vi(X?)dBi 
Jo 

are understood in the sense of Young's integration; see OS an d H2] ■ But if H > | the integrals 
that appear in (|l.lf) are understood in the rough paths sense of Lyons (see [5]). For the convenience of 
the reader, we included in an appendix at the end of this paper some results of rough paths theory that 
are used in our proofs. 

By using |5| and Theorem 6.3.1. pp. 179 of ^31> it is possible to show the existence and the uniqueness of 
a process {X^°) t >o solving Observe that from the change of variable formula the process (Xf°) t >o 

is such that for every smooth function / : R n — » R, 

f(xn = f(x ) + E f \v t f)(Xf°)dBl. 
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We denote by C b °°(R",R) the set of compactly supported smooth functions K" R. If / e C£°(R n , 
let us denote 

P t /(x )=i: (/(*?")), t>0, 

where X* is the solution of at time t. 
Our main result is the following: 

Theorem 1. Assume H > -|. There exists a family (r^J fc>Q of differential operators such that: 
(1) J/ / 6 C b °°(R n ,R) and a; £ R", i/ien /or every TV > 0, lo/ien t -> 

AT 

p t f( x ) = Y^t 2kH (^f)(-) + o(t^ H y, 



k=0 



(2) 



(3) 



i=l 



Tf = *f3(2H,2H) + m H \ E 

i,j'=l ^ i,j=l 



if ii 



, I 2 

J / J 



where 0(a, b) = J* a? " 1 ^ - xf^dx; 

(4) More generally, 1^ is a homogeneous polynomial in the V(s of degree 2k; 

(5) If H = \ or ViVj — VjVi for every 1 < i,j < d, then 



T k kl2k (E^ 2 ) 



Remark 2. The proof of this Theorem relies on the explicit bound of modulus of continuity of the ltd 
map (see the appendix). In the book of |13|. this is done only for sample paths with p finite variation, 
with 2 < p < 3,. This is the only reason why this Theorem is stated for H > ^, but it certainly also holds 
true for H S 



Remark 3. In the case of Brownian motion, some more precise results are available in pQ, [2] and 0]. 
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2. Commutative case 



In this section, we investigate the simplest case which is the commutative case. More precisely, we assume 

throughout the section that the Lie brackets [Vi, Vj] — ViVj — VjVi = 0, 1 < i, j < d. 

For i = l, d, let us denote by (e* Vi )t G K the (deterministic) flow associated with the ordinary differential 

equation 

§ -«<«>■ 

Proposition 4. The flow <£> t associated with equation M.l)) is given by the formula 

$ t = e^ o...oe^< 

Proof. Observe first that since the vector fields Vi 's are commuting, the flows (e tVi )tgR are also commuting. 
We set now for (x, y) G M™ x M d , 

F(x,y) = (e yiVl o ...oe VdVd ) (x). 

By applying the change of variable formula, we easily see that the process (e B * Vd Xo ) is solution of 

V J t>o 

the equation 

d(e B ? v «(x )) = V d (e B ? v *(xo))dBf. 
A new application of Ito's formula shows now that, since Vd and Vd-\ are commuting, 

d(e B ?~ lv «-i(e B ? v «x )) = 

V d ^ (f^'-^e^xo)) dBt 1 + V d {e Bd ^y^{e B ^x Q )) dB?. 
We deduce hence, by an iterative application of the change of variable formula that the process (F(xo, B t ))t>o 
satisfies 

d 

dF{x ,B t ) = Vi(F(x Q ,B t ))dBi. 

i=l 

Thus, by pathwise uniqueness for the equation l]l.ip . we conclude that 

X*° =F(x ,B t ), t>0. 

□ 

Remark 5. Observe that the expression 



e v ^ o • • • o e v * B ? 
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is actually defined for every H £E (0,1). Therefore, in the commutative case, it makes sense to define 
solutions of stochastic differential equations driven by fractional Brownian motions without restriction on 
the values of the Hurst parameter H , and without using rough paths theory . For instance, solutions of 
one- dimensional equations like 

dX t = o-{Xt)dB t 

are well defined for any value of the Hurst parameter, see for instance |14j . 
Corollary 6. For any smooth f : R™ — > K, 



E(/(Xf)) = (cxp fy H pvA /) (.to) 



That is, the function 

<p(t,x)=E(f(X?)), 

satisfies the partial differential equation 

i=l 

associated with the initial condition 

<p{0,x) = f(x). 

Proof. Observe first that from Ito's formula for the fractional Brownian motion, see [7] 
E^e v ^f(xo))^f(x ) + Hj\ 2H - 1 E^ v ^V l 2 f(xo)) ds. 

Therefore, 

E(e^/(^o)) = (expQ^ 2 ) /) (x ). 
It has been seen that (Proposition 0} 

$ t = e^o...oe^. 

Thus 



E(/tX? )) = fexp (^t 2H fy^ fj (x ). 



□ 



Remark 7. So, in the commutative case, there is a Feynman-Kac type formula for solutions of equations 
driven by fractional Brownian motions. It shall be shown later that this type of formula only holds in 
the commutative case. 
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Example 8. Let us consider a one- dimensional stochastic differential equation of the type 

(2.1) X?°=x + f <r(X*°)dB s 

Jo 

where a : M — > R is a C°° bounded function and B is a fractional Brownian motion with Hurst parameter 
H G (0, 1). Then, the function 

satisfies the partial differential equation 

d V ff ,2H-l 2/ 
-Q- t = Ht 

associated with the initial condition 

lf (0,x)=f(x). 

3. Asymptotic development in small times of P t 

We now study the generic case of non-commuting vector fields. Throughout this section we assume H > | 
and introduce the following notations: 

(1) 

A k [0,t] ={(*!,..., t fc ) G [0,t] k ,h < ...<t fc }; 
(2) If I = (ii, ...i k ) G {1, d} fc is a word with length k, 

[ dB 1 = [ dB\\...dB\l. 

JA k [0,t] J0<t 1 <...<t k <t 

Theorem 9. For f G C£°(R n , M) and TV > 7 when t -> 0, 

= /(*)+ X;* 9 ** E pv-^/h*) / ds / +o(^+ i ^), 

fc=l /=(u,..^) • /A2fc [°' 1 l 

and 



fc=i /=(n,...i 2 ,) V-Wi] y 



)• 



= (<l,...<2fc) 

Proof. Let us denote by -B m the sequel of linear interpolations of B along the dyadic subdivision of mesh 
to, that is if = i2- ,n for i = 0, 2 m then for t G [i™, i^], 

i? t m = B t 7 + 2 m (t - C) - B£) . 

Consider now the equation 

(8.1) ,Y t m ' x = x + J2 f V t (Xr' x )dB^ m . 
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The process X m ' x , defined in (|^ . 1|) has Lipschitz continuous sample paths. Let p > -jr. According to 
Theorem 5 of 9 , (X T t n ' x , t € [0,1]) converges to (Xf, t S [0,1]) in the distance of p variation (see 
Appendix (|6.5[) for the definition of this distance) . 

Let / be in C£°(R™,R). Using 2N + 1 times the change of variable formula 



f(Xm = f(x)+J2 I (Vrf)(X™>*)dB 



we obtain 

2N+1 



(3.2) f(Xr>*)=f(x)+ £ E (Vn...ViJ)(x) [ dB 1 



I=(n,-,*2N + 2) Jo<ul< - <U2N+2<t 
By taking the expectation we obtain therefore: 

2N+1 / \ 

(3.3) E(f(XD) = f(x)+Y / E (Vn-VUMEi dB 1 '™) 

+ E E (/ W 1 .»^ +s /)W/)<iB*\...d^+» ] • 

/=(il,...,i2N+2) \ J 0<n 1 <...<t l2JV+ 2<t / 

Since / is continuous and bounded the left member of (US. -ill converges to ~P t f{x) when m goes to infinity : 

(3.4) lim E(f(X^ x )) = Pt/(z). 

m — >oo 

Let now B m = (LB™' 1 , ...,M m - 2N+1 ) be the smooth functional over B m in the sense of Definition 3.1.1 
page 30 of ^Sj (see also the Example 1211 of the Appendix EJ. For k < 2N + 1, I = (ix, ...,«&) we have 
therefore 

dB m ' 7 = B m < 7 . 



/A fc [0,t] 

According to Theorem 4 of [H], B" 1,2 converges, in the distance of p variation p > jj, almost surely and 
in L 2 to the geometric functional denoted by B 2 . For k < 2, we thus have 



L 



dB 1 = B 2 ' 7 . 



/A fc [0,t] 

According to Theorem 3.1.3 of or Theorem [25] of the Appendix El the geometric functional B 2 has 
a unique extension in C , P (T ( - 2N+1 \W 1 )) denoted by B 2W+1 , and for k < 2N + 1, I = (h, ...,i k ) almost 
surely 

dB I =B 2N+lJ = Um [ dB rnJ 
A k [0,t] m ^°° JA k [0,t] 



8 FABRICE BAUDOIN, LAURE COUTIN 

Since J A h< t i dB m I belongs to the feth Wiener chaos of B, the limit also holds in L 1 according to 0. 
That is, for k < 2N + 1, I = (h, ...,i k ) 

Elf dB 1 ) = E (M 2N + 1 ' 1 ) = E | lim / dB mI ) . 

\^A*[0,t] / \ m ^°° J A"[0,t] J 

Let us now observe that by symmetry for k an odd integer 

E ( / dB m I J = 0. 

\JA k [0,t] J 

Using scaling property of fractional Brownian motion, we have for k an even integer 

E ( f dB 1 ) = t Hk E ( f dB 1 ) . 

\JA*[0,t] J WA fc [0,l] / 



Define for k = 1,...,N 



(3.5) r 



= lim Y El f dB m A V h ....V i2k . 
Then, the first sum in the right member of (|3.3|) converges to Ylk=i t 2hH Q?Ff)(%)i that is: 

2N+1 / \ AT 

(3.6) Hm £ E (^-^/)WE / dB I < m )=Y J t 2kH {T%f){x). 

According to Theorem 17 of |Hj , for any q > 1 and p > -jj there exists a random variable C p belonging 
to L q such that for any m, k = 1,2, I = and (s, t) G A 2 1 j, 

\B#\<C p \t-s\ k '>. 

In what follows and k may vary from lines to lines. According to Theorem 1371 there exist 9 > 1, k > 1 
( depending on a;, V^, i = and p), such that the geometric functional Z" 1,2 over (B m , X m ' x ) is 

controled by u>(t, s) = k(1 + C®)\t — s\ for any m G N. Then, according to Theorem |3H applied to 

o(6, «)(«,«;) = (v,.f(x)v) V(b,x) G R d+n , V(u,iu) G M d+ " 

there exist > 1, k > 1 ( depending on x, f, Vi, i = 1, ...,d and p), such that the geometric functional 
Qm,2 over ^ j^x m - x )) is controled by uj{t, s) — k(1 + C®)\t — s\, for any meN. Therefore, according 
to Theorem 1281 there exist 9 > 1, k > 1 ( depending on x, /, K:, i = 1, ...d, N and p), such that for any 
/ = 12N+2), rn 6 N, and i 



(3.7) 



0<«l<...<tl2JV + 2<* 



< k(1 + C p 9 )i( 2JV+2 '/P. 
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By taking the expectation of each members of (|3.7|) and using the fact that C® belongs to L 1 , we 
deduce that there exists a constant k depending only on x, /, Vi, i = 1, ...d, N and p, such that for 
I = (ix, 12N+2), m 6 N, and i 

(3.8) E(| / f{X^)dB^K.dB^^S<^ N+2),P - 

J 0<Ui<...<U2N + 2 

Finally, by taking the limit of each term of 1)3. 3|l and by using (|3.4|) . (|3.ti[) and (|3.8[) we get Theorem 

El □ 

4. Expectation of iterated integrals of the fractional Brownian motion 
In this section we analyse the coefficients appearing in the differential operators 

1 1, — 



We shall need the following lemma 



Lemma 10. Let G = (G±, G2fc) be a centered Gaussian vector. We have 



E(G 1 ...G 2k ) = ] ^ n E ( G ^)G.( 2i -i)), 
<je6 2t i=i 

where &2k is the group of the permutations of the set {1, 2k}. 
Theorem 11. Assume H > g. Let L = (ix, ...,«2fe) be a word, then 



Ell dB 1 



A 2fc [0,l] 



21, 



iL(2H-l) fe £ f n^(»-D.v(»o 1^(20-^-1) \ 2H - 2 d tl ...dt 
where 5i_j is the Kronecker's symbol. 

Proof. Let us first hint how this formula works with a heuristic argument. We have 

E ( / dB 1 ) = f E (dB 1 ) . 

\JA»[0,1] / JA 2k [0,l] 

But, from the Lemma ITUl 

1 k 

E (dB 1 ) — — — j: V nEte M2i, dR M2i - 1) V 
By using the covariance function of the fractional Brownian motion, we get therefore 
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which leads to the expected result. 

We now turn to the rigorous proof. Let us again denote by B m the sequel of linear interpolation of B 
along the dyadic subdivision of mesh m and recall that 

Elf dB 1 | = lim E ( / dB" 1 ' 7 | . 

\JA™[0,t] J " wo ° \JA^[0,t] J 

^From the Lemma llOl we have 

/ dB m,i ) = J_y TTe ( — ——j — — ) dh...dt 

J A 2„ [0>t] J k\2 k ^J A2kM ±± y dt a(2l) dfyai-i) J 

If t a (2i) e [tT^T+il and *«7(2/-i) G [t] n ,t" l +1 [, using the expression of B m , we get 



2k- 



<&-<;(%) dt a (2l-l) 



E ( — 31 1 = ^-^o 2 > fOT 4 = Ji 

= ^ (21 -i„M 2 o 2miE ( A rSA™ 1J B), for j = i-l; 

= <W 1)>i<r(2i ,if(2H - l)2 2m / / |ar - y^" 2 ^ for j < f - 1. 

Here, A™B is the increment of B between if and that is, A™B = B{tf +1 ) - B(t™). By using 
Cauchy-Schwartz inequality, we have 

|E(Af 5Ai_i5 ro )| < 2~ 2mH 
and the result follows from the Lebesgue convergence dominated theorem. 

□ 



Lemma 12. Let H > i. 



2 ' 

2H-2- 1 



(fa -tiY^-'dtidti = 



A 2 [o,i] 2H(2H-l) 

(*4 - <3) 2H ~ 2 (<2 - h^^dhdhdhdh - ' { 2IL ' 1U ' 



A4[0il] ^ ^ —™ 4ff(2H-l)2 

(t 4 - h) 2H - 2 {t 3 - t 2 ) 2H - 2 dt 1 dt 2 dt 3 dt A ! 



. ... . 2H-1U + \2H-2j+ j+ j+ j+ _ 1 PC*H,2H) 

A 4 [0,1] 



(f 4 - fcr-'fo - ^r-dtx*^ - 4g (4 g -i)( 2g -i )2 - 4g (2 g -i )2 



Proof. The result follows by direct but tedious computations. 
First of all, we recall that the beta function 



l 

/3( a ,b) = I x a -\l- x) h - l dx 

satisfies: 



(3(a + l,b) = |/3(o,6 + l) 
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0(a,b) = 0(a + l,b) + 0(a,b + l). 

(1) For the first integral: 



(2) For the second integral: 



/ (*4 - t 3 ) 2H - 2 (t 2 - t 1 f H - 2 dt l dt 2 dt i dU 
Ja 4 [o,i] 

/>1 />t4 /-t3 />t2 

/ / / / (t4-t 3 ) 2H ' 2 (t2-t 1 ) 2H - 2 dt 1 dt 2 dt 3 dt4 
Jo Jo Jo Jo 

0(2H, 2H) 



4H(2H-1) 2 ' 
(3) For the third integral 



'A 4 [0,1] 
1 

~2H - 1 
On one hand, 



/ {U - h) 2H ~ 2 {t 3 - t 2 ) 2H ~ 2 dt 1 dt 2 dt 3 dU 

JA 4 f0.11 



'' t3 {tf- 1 (U h) 2 "- 1 ) (t 3 - i 2 ) 2ff - 2 dt 2 dt 3 dt 4 
o Jo Jo 



X — f £ J** tl H -\t 3 t 2 ) 2H - 2 dt 2 dt 3 dt A 



2H 

1 



{2H-lfJ Jo 
1 

~8H 2 (2H - l) 2 
On the other hand, 

-1 rU r t 3 



1 /-t 4 

£4 ^3 dt 3 dt/± 



10 Jo Jo 

fl /•! /•! 



(*4 - t 2 ) 2H -\t 3 - t 2 ) 2H - 2 dUdt 3 dt 2 
U, r.2) 2H - 2 ((l-t 2 ) 2H -(t 3 -t 2 ) 2H )dt 3 dt 2 

'dt 2 



2H - 1 
1 

= 2F- 1 
1 



/0 Jt2 ^t 3 

f 1 /•! 



2#(2# - 1) 7 

= 2ff(2ii-l)2 / (1 ^-^2 - 2J y (2J y_ 1 1)(4ff _ 1) / (! - ^) 4if -^ 

By putting things together, we obtain the expected result. 
(4) For the last integral, the computation follows the same lines so that we do not enter into details. 
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□ 



As an immediate corollary, we deduce 
Corollary 13. Assume H > i. 



i=l 



i,j=l v ' i,j=l \ \ ) / jj—! 

It is interesting to observe that the previous corollary makes sense and actually also holds true for H > j. 
Theorem 14. Assume H > 4, i/ien i/ie conclusions of Corollarv \l?-A are still true. 

We devote now the end of the section to the proof of this theorem. The trick is to perform transformations 



dB mI 

A 2k [0,t] J 

before passing to the limit. 
Showing that for H > j, 

i=l 

is easy. Indeed, let (ii,h) £ {l,---,^} 2 such that i\ ^ «2- It is easily seen that B m ' 11 and B m t2 are 
independent, thus 



E | / dB m I = 0. 

A 2 [0,1] / 



Now, if i e {1, d} and / = (i, i) then 

E 



(7 dB m A = l -E({(B m \l)f) = \. 
yA 2 [o,i] / 



Since rf^ = J2i=(i j) lim m ^oo E y J A 2[ ^ dB mJ J VjVi : we get the expected result. 

Let now / = (i, j, k, I) G {1, c?} 4 . Recall that -B m is absolutely continuous with respect to the Lebesgue 

measure, and the fourth iterated integral is 

J r m f dB m >\ dB m >i dB m > k dB m \ 

dB ' = — — [ui) — (U2) — ("3j— ; — {UA)duidu2du^du4 L . 

A 4 [o,i] yo<«i<«2<«3<«4<l ""u du du du 
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Applying Lemma ITOl yields 
E 

where 



[ / dB 1 '" 1 ) = 5 hl S Lk A mA + S t . k 5 u A m > 2 

V-Mo.i] J 



m,3 



JO 



du du 



-A ' = / du 4 / du 3 E(-r— faa)^— M) / dw 2 / d Ul E(— — (u 2 )— — (ui)); 



(i 



A = du 4 du 2 E(-—(u 2 )-—{u 4 )) / du 3 / d Ul E ^— u 3 )-^ Ul )); 
Jo Jo du du J U2 J 



du du 



Ui , „, dB m , , dB r ' 



A mi = du 4 du 1 E{——{u 4 )——(u 1 )) / du 3 / du 2 E{ — — (ua)— — (u 2 )). 

du du J Ul J Ul du du 



'0 Jo 
Therefore, we have to prove that 

(4.1) 

(4.2) 
and 
(4.3) 

The proof relies on several Lemma 



U4 pU$ 



,dB m , dB r 



lim A mA = —(3(2H,2H); 

■n—>-oo 4 



lim A m ' 2 = — — — — r- - - 
m^oo A(AH - 1) 4 



/3(2if, 2#) 



l im A" 1 ' 3 = 2H \ 
m^oo 8(4H - 1) 



Lemma 15. For any a < H, the random variable 
sup \\B m — B\\ a = sup sup 

m meN(s,()6[0,l] 2 ,s<l 

/las a finite exponential moment. 



\{B(t)-B m (t)]-{B( S )-B m ( S ) 
\t-s\" 



Proof. According to Kolmogorov lemma, the random variable sup( s t)e[o i] 2 s<t " B< j^_T s ^' S ^ is finite. Since 
B is Gaussian process, according to Theorem 1.2.3 of |12| , it has a finite exponential moment. Then the 
Lemma is a consequence of results on linear interpolation of Holder functions. □ 



Lemma 16. For H > j, 



H 



lim A m ' L = — p(2H, 2H) 

ra — >oo 4 



Proof. Integrating with respect to Ui and u 2 in the expression of A" 1 ' 1 leads to 



A™' 1 = - / du. / du 3 



du 3 du 4 



E(B m {u 3 )B m (u i ))E(B m (u 3 ) 2 ). 



14 



FABRICE BAUDOIN, LAURE COUTIN 



Using Fubini's theorem, we integrate first with respect to 1*4: 



A" 1 ' 1 = i jf du 3 E([B m (l) - B m (u 3 )]^ 1 B m (u 3 ))E(B m (u 3 ) 2 ); 



du 3 ^-E(B m {u 3 )B m {l))E(B m {u 3 ) 2 ) - \E(B m (l) 2 ) 2 
du 3 4 



Now, using the expression of B m we obtain that for u 6 [i™>^i+i[' 



-^—E(B m (u)B m (l)) = 2 m E(A™BB(f)) 
du 



2H 



2 

H2 



i\t? + ,r-\i-tr + ir-(\t?\ 

^W 2 "- 1 + \l-r\ 2H - 1 ]dr. 



|2ff \ 



Then, using Fubini's theorem, we get 



where 



duE(B m (l)-^-B m (u))E(B m (u) 2 ) = H / Ir 21 *' 1 + |1 - r\ 2H - x ]a m {r)dr 
du Jo 



2"' -1 



a m (r) = 2 m / E(i? m ( U ) 2 )d U l [tm , tm+i[ (r). 

For all r £ [0, 1], a m (r) converges to E(B(r) 2 ) = r 2H . Following Lemma IT51 it is bounded uniformly 
m and r. Then, using dominated convergence Lebesgue theorem we obtain 



f 1 d 
lim / duE(B m (l)—B m (u))E(B 

i^+oc J Q du 



m {u) 2 ) = H I y"- 1 + |1 - r\ ztl - l ]r 2tl dr 
Jo 

I + f/3(2if,2F) 



because /3(2#, 2ff + 1) = |/3(2F, 2#). We conclude that 



lim A" 1,1 = —j3(2H,2H). 



Lemma 17. For H > j, 



lim A m ' 3 = 



2# - f 
8(AH - 1) 



Proof. Recall that 



A m - 3 = / du A I duiE(^(tt4)^(«i)) /""dug r du 2 E(^(u 3 )^(u 2 )). 



du 



da 



du 



du 
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Integrating with respect to ui and U3 we obtain 

A m > 3 = -j du 4 J d Ul E( — M^-( Ul ))E([i? m M - i? m (ui)] 2 ). 
We introduce the indices i,j such that 114 € [^'^i+iL an< ^ w i ^ Recall that 

= 2 2m E(A™BA™!B) for j = i - 1; 

= H(2H - l)2 2m // |a; - y\ 2H - 2 dxdy for j < j - 1. 



Then, we split A m ' 3 into three parts 



where 



2"' -1 



A 3,i, m= J2 2 4m - 2Hm / / dv{u-v) 2l &{&™B 2 )- 

2 '"- 1 rtf+x ft? 

A 3,2, m = 2 2m E(A^SA^ X B) / du / d«E([B m (u) - £ m (v)] 2 ); 



and 



A 3 > 3 ' ro = If (2tf - 1) V J2^ 2m [ t,+1 dx f 3+1 dy{x - y) 2H ' 2 f' +1 dudvE{[B rn (u) - B m {v)) 2 ). 



i=2 j=0 

We have to prove that for k = 1,2 

(4.4) lim A 3 ^" 1 = 

rn— >oo 

and 



(4.5) lim - 2,1 ' 1 



4(4# - 1) 

First we prove 14.411 . Integrating with respect to u and v yields 

^3,1, »™ _ J_ 2~ 4Hm 

~~ 12 ^ 



i=0 

Since H > j, we deduce that lim m _ 00 A 3 ' 1 '" 1 = 0. 
According to Lemma ll 51 we have for a < H 

l^ 1 '" 1 ! < , 9 +Z 9 -^ 2( 1 -^- 2 ")^E( S up snp |J3 y~^ )|2 ). 
(2a + l)(2a + 2) m ( s , t )e[o,i] 2 ,s<t I* - T° 
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Since H > |, we deduce that linim^oo A 3 ' 2 '" 1 = 0. 

Now, we prove (|4.5|l . Indeed using Fubini's theorem, we have 

A 3 ' 3 ' m = H{2H -1) ( dx [ (x-y) 2H - 2 a 3 > 3 ' m (x,y)dy, 
Jo Jo 

where for < j < i - 1, and (x,y) e [C^i+iM*™.^! 

a 3 > 3 ' m (x, y) = 2 2m f * +± du f ]+% dvE((B m (v) - B m {u)) 2 ) 



and a 3 ' 3,m (x, y) = elsewhere. For all x,y a 2 ' 3 ' m (x,y) converges almost surely to E([B(x) — B(y)] 2 ). 
Moreover, since \f^_ 1 — i™| < 2\x — y\, a 3 ' 3 m is bounded by 

la 3 ' 3 '™^)! < 2E(supsup lB{ ^f} t)l ) 2 )\2(x - y)\ 2 « . 

m s<t |£ S\ 

Using Lemma^|and dominated Lebesgue convergence theorem, we take the limit when m goes to infinity 
for H > \: 

lim A 3,3 '" 1 = H(2H -I) [ dx [ (x - y) 2H - 2 E{[B(x) - B{y)] 2 )dy 

rn-nx) J J 

2H - 1 



4(4ff - 1) 



□ 



Finally, 

Lemma 18. For H > \, 

nl~'oo" 4(4#-l) 4 

Proof. We know that if the vector fields commute then 



lim A m ' 2 = — — - —/3{2H,2H) 



2 



Therefore 

lim A m ' 1 +A m ' 2 +A m ' 3 = -. 

m— H-oo 8 



□ 
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5. Application to the study of invariant measures 

Consider the stochastic differential equation on R™ 

(5.1) dX t = Vi{X?°)dBi 

where the Vi's are C°°-bounded vector fields on M. n and B is a d dimensional fractional Brownian motion 
with Hurst parameter H > g. 

Proposition 19. Assume that [i is a probability measure on K™ that is invariant by 1)) . In the sense 
of distributions, we have for every k > 1, 

(rf)*M = o, 

where (T^)* is the formal adjoint ofrj?. 

Proof. Let / e C b °° (W l , K) . We have for every N > 0, when t -> 

E ( f f(X?Mdx)) = jh t 2kH [ (T«f)(xMdx) + o(t^»). 
But since fi is invariant, 

E (7 fiXfMdx)) = [ f(x)^dx). 
The result follows therefore. □ 

6. Appendix : On rough path theory 

We recall some definitions and results on rough path theory, see [EJ. Indeed, we precise how all the 
constants appearing in the continuity Theorem of the Ito map depend on the vectors fields and the 
control of the paths. 

6.1. Basic definitions and properties. We work on V = M. d endowed with the Euclidean norm. The 
tensor product is V® k = V ® ... ® V (of k copies of V) endowed with a norm \.\k compatible with the 
tensor product that is 

\t®ri\k+i < \S\k\v\i, V^V® k , V V eV® 1 . 
For each n^N, the truncated tensor algebra T^ n \V) is 

n 

T {n) {V) = ^2 v ® k ' v ®° = K - 

fc=0 
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Its multiplication is 

k 

(£®r,) k =J2^ ®V k ~ j , k = 0,...n, V£, i)6r w (F). 

3=0 

The norm |.| on T^(V) is denned by 

n 

The pair (T (n) (V), |.|) is a tensor algebra with identity element (1, 0, 0) and for £, 77 6 T (n) (V), |^<8>7?| < 

l£IM- 

The tensor algebra T^(V) is 

OO 

r (oo) (^) = £ ^® fe , y®° = m. 

We use Ajp jij to denote the simplex {(s,i), < s < t < T}. Recall that a control w is a continuous, 
supper additive function on A^ T j with values in [0, +oo[ such that u>(t, t) — 0. Therefore 

Oj(s, t) + 0j(t, u) < U}(s, U) V(S, t), (t, u) 6 Ar Q)T j. 

Definition 20. j4 continuous map X from the simplex A? T -, into a truncated tensor algebra TW(V), 
and written as X s> t — (X^ t , .., X^ t ) with X k t e V® k for any (s,t) £ A^ 0T j, fc = 1, ...,n is caZ/ed a 
multiplicative functional of degree n (n £ N*) if 

(6.1) X s . t ®X t , u = X s , u , V(s,t), (t,u)eAf 0T] , 

where the tensor product <g> is taken in T^(V). 

Equality 1)6. lfl is called Chen identity, although it appears long before Chen's fundamental works in which 
a connection is made from iterated path integrals along smooth paths to a class of differential forms on 
a space of loops on manifold. 

Example 21. Let x : [0,T] — > V be a continuous path. If x is Lipschitz path, then we may build a 
sequence of iterated integral X k t — J s<tl< <t fc <t^ x ti ® ■■■®dxt k . In this case identity f \b'.l}) is equivalent 
to the additive property of iterated path integrals over differents domains. 
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Definition 22. Let p > 1 be a constant. We say that a map X : [0, T] — * T^ n '(V) possesses finite p 
variation of 

\Xy Vi = l,...,n, V(s,t)eAf 0>T] 

/or some control to. 

Definition 23. A multiplicative functional with finite p variation on T^ P ^(V) is called a rough path (of 
roughness p). We say that a rough path (of roughness p) X in T^''(V) is controled by uj if 

|X; t | <w(«,t) </p , Vi = l,..,b], andV(s,t)GAf 0iT] . 
The set of all rough path with roughness p in 

T(W)(y) w m be denoted by fl p (V). 

Definition 24. A smooth rough path X is an element of fl p (V) such that there exists a Lipschitz 
x : [0, T] — -> V Lipschitz such that 

X k s t = I dx H ®...®dx tk , k = l,...,\p\, V(«,t) € Af 0iT] . 

J S<tl<...<tk<t 

6.2. Extension of rough path. The following theorem shows that the higher (than \p\) order terms 
X k (k > [p]) are determined uniquely by X 1 (i < [p]) among all possible extensions to a multiplicative 
functional wich posses finite p— variation. 

Theorem 25. Let p > 1, and let X : A^ T , — > T( n )(V) be a multiplicative functional with finite p 
variation so that 

|Xj it | <aj(s,t) l/p , Vi = l,..,n and V(s,t) G Af T] , 

/or some control ui. If n > [p], £/ien we may uniquely extend X to be a multiplicative functional in T^°°\V) 
with finite p variation. Moreover, if to is a control such that 

(6.2) \XIA< W( ^ P , V» = l,...,[p], andV(s,t)e Af 0iT] , 

where [3 is a constant such that 

r=3 

i/ien |6'.ii|) remains true for all i > [p] . 

The extension of a rough path X to a higher-order multiplicative functional is continuous in p— variation 
distance. 
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Theorem 26. Theorem 3.1.3 v 39 of [H] 

Let X and Y be two rough path of roughness p and let (3 be a constant such that f3 > 2p 2 [l+X^3(;rz2) p ]• 



If uj is a control such that 



(6-3) \Xl t \, |y;,tl< ^y , Vi=l,-,[p], andV( S! t)GAf 0iT] , 



(6-4) Kt-^tl <g ^|y^ ; V» = l,...,[p], andV(»,t)6Af 0i r,, 

</ien |6'.ffi and jo'.^| ) TioZd for all 



i. 



6.2.1. Almost rough path. In this section we give a method of constructing rough paths. 

Definition 27. Let p > 1 be a constant. A function X : A 2 — > r^ p ^(V) is called an almost rough path 
( of roughness p) if it is of finite p— variation, Xg t = 1 and for some control uj and some constant 6 > 1, 

\(x s , t ®x t , u y -x StU \<oj( s ,u) e 

for all (s, t), (t, u) G A? Q T j and i = 1, [p]. 

The following theorem justifies the name of almost rough path. 

Theorem 28. If X : A 2 — > T(W)(y) is an almost rough path of roughness p, controlled by uj and 6 then 
there exists an unique rough path X (with roughness p) in T^\V) such that 

\X\, t -X'l t \ < Kiuj(s,t) e , VI < i < [p], V(a,i) G Af 0T] , 

/or i*Q defined by induction 

Kq = maxuj V 1, 
00 / a 



r=3 



r-2) ' 

1 + J2 {?K l J p K k+1 ^ + KiKk+x^Kg 



i=i 



The following theorem shows that in fact the map X — > X is continuous. 
Theorem 29. Let X andY be two almost rough paths of roughness p in 

T(M)(V), of which controlled 

by a control uj, that is 

\Xi t \, \Yl t \<uj(s,ty/P, Vi = l,...,\p], V(s,t)G Af 0iT] . 
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and for some 6 > 1 \(X s>t (g) X t:U ) 4 / p — X\{u < w(s, u) e /or aZZ (s, i), (t, u) e Ajp T j, i = 1, [p], wiZZi i/ie 
same inequality also holding for Y. Suppose that 

\Xi t -Yi t \<eu(s,ty/P, Vi = l,...,[p], V(M)eAf 0i31) 

l^, t -^l<5i(^Kt) i/, ) Vi = l,...,[p], V(s,i)eAf 0;T] ; 
where the Bi are defined inductively by 



. r=3 



r-2 

B k+l {e) = e + 2,^A k {r 1 e) 



i/p 



r-2 



e A if e , 



r=2 



and 



i4 fc (r, e) - min{if e ( — -) e (l + 2 V if fe+1 _,if* /p + i^iffc+^if*), 



( )(fc+i)/ P[ £; [B . (e ) ( i + ^ fc+1 _ 2 ) + # e ] + e]}) 



i=l 



r-2 

if = max if j 

»=i,-b] 



6.2.2. Spaces of rough path. Let C(Aj? T j, T( n )(V)) denote the set of all continuous functions from the 
simplex Ajj, into the truncated tensor algebra (T(")(V), |.|). If X G C(Aj? ^ T( n )(V)), then we may 
write 

X Stt = (X° Sit ,...,X^ t ), V(M)eAf 0;T] , 

where t S V® 1 is the ith component of X (also called the zth level path of X). The subset of the 
functions of C(Af Q T] , T^(V)) such that X° s t = 1 is denoted C {Af Q T] , T^(V)). 

Definition 30. A function X G C (Aj* T ^,T^ n \V)) is said to have finite total p-variation if 

supVl-Xj t !"/*<(», i = l,...,n, 
w/iere supc r«ns ewer aZZ /Zraie subdivisions of [0,T]. 
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Proposition 31. Let p > 1 be a constant, and let X 6 Co (A? r i , T^ 1 ' (V)) satisfy Chen's identity \6.1]) 
( i.e. X is a multiplicative functional inT^(V) of order n). If X has a finite p-variation, then 

n 

w(s,t) = E SU P E \Xu-uti\ p,i , V(*,t) G Af 0)T] 
»=i i 

is a control function, and 

\Xl t \<cj( S ,ty/P, Vi = l,...,n, V(«,t)GAf 0i31 . 

Let Co,p(A? T] ,T^(V)) denote the subspace of all X E C (A^ Q T] ,T {n \V)) with finite p variation. It is 
clear that Cb, p (A? T ^T^ n \V)) is a metric space. The p variation metric d p on Co, p (A^ T ^T^ P ^(V))) is 
defined by 

(6-5) d p (X,Y) = max sup ( E I^Li.t, " K-ut 

t=l,...,[p] D ^ 

The space (f2 p (V),d p ) is a complete metric space. 

Definition 32. Geometric rough paths with roughness p are the rough paths in the closure of smooth 
rough path under the p— variation distance. 

The space of all geometric rough paths with roughness p is denoted by Gfl p (V). 

6.3. Integration theory degree 2. Let p e [2,3[. Let W be R™, n > 1. Let a : V — > L(V,W) be a 
function which sends elements of V linearly to W- valued one-forms on V. Suppose that a posseses all fcth 
continuous derivatives d k a up to the degree 3 and denote a 1 = d l a, i — 1,2, 3. 

Let X G Q, p (V) and let X s t = (1, X^ t , X^ t ). The almost rough path which defines the path integral 
/ a(X)dX is Y G C (A2 0T] ,T( 2 )(W^)) where Y s>t = {l,Y s ) t ,Y* t ) and 

Y lt = al i x o.s)- x lt +a 2 ( x a, s )- x s,t> 
Yl^a^XlJ^a^XlJ.Xlt. 

Theorem 33. Theorem 5.2.1 and remark 5.3.1 of [T3] 

Let a : V — > L(V,W / ). Suppose that a possesses all kth continuous derivatives d k a up to the degree 3 
and 

\d l a(0\ U vx...xv.,w)<M(± + \t\), i = l,...,3, V£ G V. 




SDES DRIVEN BY FBM 23 

Assume that X S fl p (V) is controlled by u>, namely 

\Xi tt \<u(8,ty^, i = 1,2, V(«,t)e Af 0iT] . 

T/ien y is an almost rough path with roughness p in T^ 2 \W) with control lo and 9 = 3/p i.e. there exists 
a universal constant C such that 

\(Y s , t ® Y t . u ) 1 - Y s \ u \ < CMlu(s,u) 3 /p. 

Definition 34. Let X e Q p (V). Then the integral of the one-form a against the rough path X, denoted 
by j a{X)dX, is the unique rough path with roughness p in T^ 2 \W) associated to the almost rough path 
Y e C (Af 0iT] ,T( 2 )(WO) where Y„ tt = (l,Y^ t ,Y^ t ) and 

Y s \ = aHxl s ).Xl t+a 2 (Xl s ).Xl t , 
Yl t = a\Xl s )®a\Xl s ).Xl t . 

Theorem 35. Let a : V — > L(V,W). Suppose that a possesses all kth continuous derivatives d k a up to 
the degree 3 and 

|d i a(OlL(vx...xv,wo<M(l + |^|), i = l,...,3, V£ e V. 
Assume that I, le & P (V) is controlled by u, namely 

\Xl it \, \Xi it \<w(8,ty'>, 2 = 1,2, V(s,i)eAf 0;T] , 

and 

\Xi it -Xi tt \<eu(8,tf*, i = 1,2, V(s,i)eAf 0iT] . 

T/ien 

f a{X 0tU )dXl- f a{X 0tU )dXi\ < KeMu{s,uf", 

J s J s 

for all (s, t) 6 A^ T j and i = 1,2, w/iere K is a constant which is polynomial in M, maxw. 
6.4. Ito maps : rough path with 2 < p < 3. 
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6.4.1. Framework. Let V = R d and W = E n . Let / : V -» L(V r , be a function, which can be viewed 
as a map sending vector of V linearly to a vector field on W. Consider the following differential equation 
( initial value problem) 

dY t = f{Y t )dX u 

(6.6) Y = y . 

Since the integral J f(Y)dX for rough path X, Y make no sense generally we are not able to iterate 
differential equation eqn 1)6. 6[1 to obtain the unique solution directely. To overcome this difficulty, the 
idea is to combine X and Y together as a new path. We view equation 1)6. 6JI as 

dX t = dX t , 

(6.7) dY t = f(Y t )dX t ,Y = y . 

The initial condition of X is irrelevant, and therefore we simply take Xq = 0. Define / : V © W -— > 
L(V © W; V 8 W) {V ®W is the direct sum of V and W) by 

f(x,y)(v,w) = (v,f{y + y ).v), V(x,y) e V S)W, \/(v,w) 

Then eqn l|6.7(l can be written in the following more appreciating form 

(6.8) dZ t = f(Z t )dZ t . 

Given a rough path X in V, we said that a geometric rough path Z in V © W is a solution to (|6.7() if 

n y (^) = x, 

Z= I f{Z)dZ, 

where Ily is the projector on T^\V). 

6.4.2. Existence and uniqueness results. One can summarize the results proved in pages 149 to 162 of 
1 1 3 j - when the control is u>(t, s) = C\t — s\ and the vector field, /, and its derivatives are bounded in the 
following way. 

Theorem 36. Let f E C 3 (W, L(V, W)) be a vector field and let M be a constant such that 

\d\f(£)\<M, V£eW, i = 0,l,2,3; 
\d l f(0-d l f(v)\ <M\£-T,\, veW, z = 0,l,2,3. 
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Let X be a rough path in T^ 2 \V) with roughness 2 < p < 3 controlled by u>, where for a constant C p , 
w{t,s) = C p \t-s\, (s,t) eAf 01] . 

Then, there exist some constants k, and 9 depending only on M, p such that if 

(6.9) T x = k — < i 

there exists a unique Z £ Cl p (V®W) such that Hv(Z) = X and Z satisfies the following integral equation : 

Zlt = f f(Z)dZ\ i = 1,2, V(«,t) G Af 0iTi] . 

J S 

Moreover, the following estimation holds 

(6.10) \zi tt \<(±\ t - s \y/r, i = l,2, V( S ,i)eAf 0;Ti] . 
Finally, we may extend the solution to the whole interval [0, 1]. 

Theorem 37. Let f G C 3 (W,h(V,W)) be a vector field and let M be a constant such that 

\<ff(£)\<M, V£e W, i = 0,l,2,3; 

Wf{0-d l f(n)\<M\i-nl V£, ryeVK, * = 0,1,2,3. 

Let X be a rough path in T^(V) with roughness 2 < p < 3 controlled by u, where for a constant C p , 
w{t,s) = C p \t-s\, (s,t)eA[ 0jl] . 

Then, there exists a unique Z £ Q p (V (& W) such that Yl v (Z) = X and Z satisfies the following integral 
equation : 

zi tt = f f{Z)dZ\ i = 1,2, v(«,t)eAf 0il] . 

J s 

Moreover, there exist some constants k and 9 depending only on M, p such that the following estimation 
holds 

(6.11) \Zl t \< K (l + ci)(\t-s\y/P, i = l,2, V(5,i)eAf 0il] . 

Proof. For example, we may solve the integral equation beyond 7\ by replacing the initial condition y 
by Y Tl = U w (Z) Tl . Then, the solution is defined up to the time 

S 2 = 2Ti 

Moreover, we have 

|^ 2 |<2(^)^ + | yo |. 
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By an iteration procedure, the solution is denned up to the time 1. Let N be an integer such that such 
that Nk - § > 1. Then N is bounded by a polynomial in C p . Then using the Chen rules, identity . lf> 
and several times estimation (|6.10|) we obtain 

\Zi t \<K(l + Ce)V(N)(\t-s\y/v, i = l,2, V(a, t) G Af M , 

where V(N) is polynomial in N and then in C p . □ 
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